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ABSTRACT 

In this paper we present new observations of the gravitational lens system JVAS 
B0218+357 made with a global VLBI network at a frequency of 8.4 GHz. Our maps have 
an rms noise of 30 ^Jy beam^^ and with these we have been able to image much of the ex- 
tended structure of the radio jet in both the A and B images at high resolution (~1 mas). The 
main use of these maps will be to enable us to further constrain the lens model for the purposes 
of Hq determination. We are able to identify several sub-components common to both images 
with the expected parity reversal, including one which we identify as a counter-jet. We have 
not been successful in detecting either the core of the lensing galaxy or a third image. Using 
a model of the lensing galaxy we have back-projected both of the images to the source plane 
and find that they agree well. However, there are small, but significant, differences which we 
suggest may arise from multi-path scattering in the ISM of the lensing galaxy. We also find 
an exponent of the radial mass distribution of (3 ~ 1.04, in agreement with lens modelling 
of published 15-GHz VLBI data. Polarisation maps of each image are presented which show 
that the distributions of polarisation across images A and B are different. We suggest that this 
results from Faraday rotation and associated depolarisation in the lensing galaxy. 
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1 INTRODUCTION 

Few gravitat ional lens systems are as well studied as JVAS 
B02 18+357 (IPatnaik et al.||l993|), perhaps the best example of a 



lens system for which the method of Refsdal (1964 ) can be used to 



determine the Hubble parameter, Hq. Unlike many other methods, 
this is done in a single step and involves well understood, and rela- 
tively simple, astrophysics. We point out that there remain signifi- 
cant uncertainties in Ho as determined by traditional methods (e.g. 



H(y. In most cases the main source of uncertainty is in the modelling 
of the gravitational potential responsible for the image splitting and 
distortion. A lack of observational data to constrain the lens model 
and mul tiple-galaxy lens es are both significant factors in this regard 
(see e.g. Schechter| 200C ). A list of the time d elays m easured to date 
can be found in pourbin, Saha & Schechter (2002). 



Shanks et al 



2001). For example, the use of Cepheid-calibrated 
distances to determine the distance to galaxies hosting Type la su- 
pernovae has resulted in final values of Hq that do not agree at th e 
Icr (~10 per cent) level ( |Parodi et al"||200C| ; [Freedman etak|[200l[ ). 
This is despite a large overlap between the samples of SNIa. 

Whilst the gravitational lens route to Hq can be a particularly 
"clean" one, and despite many lens systems now having measured 
time delays, almost all suffer in one way or another from effects that 
significantly increase the uncertainty in the final determination of 



With JVAS B02 18+357 on the other hand, the potential ex- 
ists to reduce the uncertainty in the lens model to a level compa- 
rable to that in the time delay; this currently stands at 3 per cent 
(r = 10.5 ± 0.4 d; Biggs et al., 1999 ). Other authors fin d a value 
for the time delay consistent with this (Cohen et al. 200C ). This re- 



duction in model uncertainty is possible due to the large number 
of observational model constraints available from multi-frequency, 
multi-resolution imaging. Both images (A and B) of the z — 0.96 
background quasar/B L Lac are easily resolved with VLB I into two 



subcomponents (e.g. Patnaik, Porcas & Browne| 
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be close to isothermal. The lens systems with the largest numbers 
of model constraints are those that contain Einstein rings, such as 
B0218+357. The resolved arcsec-scale structure of the ring probes 
the lens potential along multiple lines of sight and along all az- 
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imuthal position angles relative to the lens centre. The ring (see 
Fig. ^ has been mapped at high resolution u sing c ombined MER- 
LIN and VLA data at 5 GHz ( Biggs et al. 2001 ) and the model 
constraine d using a version of the LensClean algorithm ( Kochanek 
& Narayan 1992); preliminary results can be found in Wucknitz 



( 2001 ). Finally, the deflecting mass is concentrated in a single, iso- 
lated galaxy and so the effect of tidal shear due to nearby structure 
is very small (~1 per cent, Lehar et al. 2000). This results in a rela- 
tively uncomplicated mass model compared to other lens systems. 
A VLA map of B02 18-1-357 can be seen in Fig. §. 

The main stumbling block to date in measuring Ho with JVAS 
B0218-F357 has been that, due predominantly to the small size of 
the system (A-B separation of 334 mas), it has been difficult to 
measure the position of the z = 0.6847 lensing galaxy relative to 
the lensed images accurately using optical (HST) data. For exam- 
ple, the pos itions derived from two NICMOS observations differ 
by 46 mas (Lehar et al. 200C) and this uncertainty translates into 
a large uncertainty in Ho- With LensClean we have estimated the 
galaxy position to a theoretical accuracy of a few mas, but this re- 
lies to some extent on the lens model and so could be biased. 

In this paper we present new 8.4-GHz observations of JVAS 
B02 18-1-357 taken with a global VLBI array that were designed to 
increase the available model constraints by other means. With the 
great sensitivity of a global array we hoped we might detect: 

(i) extended structure associated with the radio jet in both A and 
B. This is detected in 5-GHz data from the First Caltech-Jodrell 



survey 

to pr oduce maps (Fig. |l|) that are superior to those published in X u 
et al.. With a resolution of ~5 mas these new maps reveal a jet that 
extends out to 30 mas from the core in each image. 

(i i) a third image, as detected in, e.g. APM 08279-F5255 (Ibata 



et al. |1999|; Mwis et al.||2002|) and MG 1 131-1-0456 ( |Chen & Hewitt 



19931 ::hen, Kochanek & Hewitt 



1995), and 



(iii) the galaxy core. Emission from the le nsing galaxy has been 
detected in several syster ns e.g. B0957-I-561 (Harvanek et al. 1997) 
and CLASS B2045-I-265 (iFassnacht et alJ|1999|). A third radio com- 



ponent is detected in PMN J1632-0 033, but it is not k nown if this 
is a lensed image or the lens galaxy (Winn et al. 2002). 



The observations have also been used to investigate the mas- 
scale polarisation structure of this system. Whilst the cores of most 
uasars/BL Lacs are polarised at a level of about 2-3 per cent 
images A and B of B0218-H357 are both 
10 per cent at frequen - 



aikia & Salter 



1988 



much more highly polari sed, reaching 



cies 8.4 GHz and above (Patnaik et al. 



1993 



Biggs et al. 1999) 



Also, the polarisation position angles of the images are not the 
same, in conflict with the expected behaviour of a lens system 
which pr eserves the polarisation position angle on the sky (D yer 
& Shaver 1992[ ). It is thought that Faraday rotation in the magne- 
toionic medium of the lensing galaxy its elf rotates the positi on an- 
gles of the images by different amounts ( Patnaik et al. 1993 ]. Thus 



the extended nature of the emission in images A and B allows the 
interstellar medium (ISM) of a high-redshift galaxy to be studied 
on parsec scales. Both images also depolarise at low frequencies. 



2 OBSERVATIONS AND DATA REDUCTION 

The goals of this experiment required a combination of high sensi- 
tivity, high resolution and excellent (u, v) coverage in order to re- 
liably image complex and faint structures. For this reason we used 
a global VLBI array including the ten antennas of the VLBA, six 



Table 1. Global VLBI array antennas and performance at 8.4 GHz. Observa- 
tions were conducted in dual circular polarisation at all antennas apart from Nt, 
On and Yb where only right circular polarisation was available. 



Antenna 


Diameter (m) 


SEFE 


VLBA 


25 


307 


VLA (Y9)| 




41 


Effelsberg (Eb) 


100 


20 


Medicina (Mc) 


32 


270 


Noto (Nt)0 


32 


770 


Onsala (On)'' 


20 


1630 


Westerbork (Wb) 


94c 


120 


Yebes (Yb)'* 


14 


3300 


Robledo (Ro) 


70 


23 


Goldstone (Go) 


70 


20 



System equivalent flux density 

Inner nine antennas of VLA phased together 

Diameter coiresponding to summed geometric area of telescopes 

Only right circular polarisation was available 



from the European VLBI Network (EVN), the VLA (Y) in phased- 
array mode and the two 70-m antennas of the Deep Space Net- 
work (DSN) at Goldstone and Robledo. Details of these antennas 
are given in Table [l| with their diameters and system equivalent 
flux densities (SEED). We chose to observe at 8.4 GHz instead of 
at 5 GHz (where the surface brightness sensitivity is better) due to 
the improved resolution and because there is less chance that the 
lensed emission would be affected by scattering; there is evidence 
that the image sizes in this system may be enlarged by scattering at 
low frequenci£s_^Biggs_£^^ly_^i^^re2ara^ 

& Patnaik ( |l996b^ for a I.7-GHz map that illustrates the large size 
increase). Scattering would possibly reduce the surface brightness 
of the third image or the galaxy core which is clearly undesirable 
when searching for such weak features. 

Due to the fact that the VLA was in its largest configuration 
('A') at the time of the observations, using the full complement of 
27 antennas as a single phased array was not possible as the syn- 
thesised beam would have been significantly smaller than the sepa- 
ration between the two images A and B. Therefore we were forced 
to use only the inner nine antennas of the array, Yg, to give a large 
enough beam. It should be noted that in order to calibrate the abso- 
lute polarisation position angle, the strategy for which is described 
later in this section, three short additional observations of 3C48 
and 3CI38 were added to the inner subarray schedule. These two 
sources were not observed by any other antenna in the global ar- 
ray. The remaining 18 antennas of the VLA observed as an outer 
subarray operating in normal interferometric mode. 

The observations took place during II November 2000 to 12 
November 2000 over the time range 20:00-07:00 UT for the Eu- 
ropean antennas and 21:00-09:00 UT for those in the US. Both 
senses of polarisation were recorded with the exception of three of 
the stations where only right circular polarisation (RCP) was avail- 
able (Nt, On and Yb). Data were recorded in four 8 MHz sub-bands 
with one-bit sampling giving a data rate of 128 Mbps. Of the 19 an- 
tennas scheduled (VLBA-l-EVN-l-DSN-l-Y), only 17 took part in the 
observations; Yebes suffered from a broken data formatter whilst 
the receiver at VLBA Owens Valley was non-operational. The var- 
ious calibration sources observed are shown in Table ^ along with 
their flux densities. The data were correlated at the VLBA correla- 
tor in Socorro producing 16 0.5-MHz channels per sub-band (IF) 
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Figure 1. CJl VLBI maps of JVAS B0218-H357 at 5 GHz. Left: B, right: A. The restoring beam is shown in the bottom-left comer of each map and has a 
FWHM of 5 X 4.2 mas at a position angle of — 1?9. Contours are plotted at multiples (—1, 1, 2, 4, 8, 16, etc) of 3cr where a is the off-source rms noise in the 
map (330 fj,]y beam~^). Both maps are plotted on the same angular scale. 



Table 2. Approximate flux densities, hour angle coverage and role of the various 
calibration sources. 



Source 



Flux density (Jy) # Hour angles Calibration role 



B0234-H285 

3C84 

0()59-^581 

0300-H470 

3C480 

3C138'" 



3.0 

19.0 

1.5 

1.3 

3.2 

2.4 



29 

8 

3 

3 

2 

1 



amplitude 

bandpass/L-R delay/D-term 
polarisation position angle 
polarisation position angle 
absolute polarisation position angle 
absolute polarisation position angle 



VLA inner antennas only 



and using an averaging time of two seconds. All polarisation pairs 
(RR, LL, RL and LR) were produced. 

Data reduction was performed using the NRAO AIPS pack- 
age. Amplitude calibration was carried out using system temper- 
atures (Taya) recorded during the observations scaled by a priori 
measurements of the antenna gains. For the DSN antennas no Tsys 
measurements were available, but values were estimated from mea- 
surements of the total power. Delay residuals were calculated using 
short segments of B0234-I-285 data and applied to all sources. This 
source was then fringe-fitted, solving for delays, rates and phases 
for the entire dataset. Following this the calibrated B0234-I-285 data 
were mapped in IMAGR and self-calibrated in CALIB. The residual 
amplitude voltage gain solutions found during the self-calibration 
were large for many antennas (up to a factor of two for VLBA Pie 
Town) and so were subsequently applied to all sources in order to 
correct the initial amplitude calibration. Next the B02 18-1-357 data 
were fringe-fitted, initially with a single component located at the 
location of the brighter A component relative to the correlation po- 
sition, but later with a more realistic model made from mapped and 
self-calibrated data. The delay difference between the left and right 
polarisations was calibrated using a single scan of 3C84 data and 
two baselines to the reference antenna (VLBA North Liberty). The 



3C84 data were also used to calibrate the bandpasses of the anten- 
nas. 

Next the various steps to calibrate the polarisation were under- 
taken. Firstly, the leakage (or D-} terms were solved for using the 
3C84 data which had first been mapped and self-calibrated; 3C84 
was assumed to be unpolarised. Secondly the absolute position an- 
gle of the polarisation was calibrated using the data from the in- 
ner subarray of the VLA (Yg) in normal interferometric modeQ As 
mentioned earlier, the bright calibration sources 3C48 and 3C138 
were included in the Yg schedule and with these it was possible 
to determine the true position angle of the source 0059-1-581. Once 
this was known it was possible to calibrate the position angle of the 
global array data by rotating the RCP— LCP phase difference un- 
til the 0059-1-581 polarisation position angle agreed with that found 
with the VLA data. 

At this point the calibrated B0218-f357 data were averaged 
in frequency to produce one channel per IF and averaged in time 
over 10 s intervals. The regions around images A and B were 
CLEANed using the IMAGR wide-field mapping task and the 



^ The antennas of the phased array also produce the usual baseline combi- 
nations as a matter of course. 
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data self-calibrated in CALIB. Many cycles of mapping and self- 
calibration were undertaken, initially only allowing the phases to 
vary, but later the amplitudes as well. It was found that the sensi- 
tivity of the Westerbork (Wb) telescope was changing by up to a 
factor of four, cycling over 20 min timescales. Due to this, and the 
fact that the polarisation leakage terms were excessively large (up 
to 45 per cent), the Wb data were flagged and took no further part 
in the analysis. Much of the previous calibration was then repeated 
to make sure that the Wb data was not biasing the calibration for 
the other antennas. Data from the Noto and Onsala antennas could 
not be included when making the Stokes Q and U maps as neither 
was able to observe both right and left circular polarisations. 

Unlike the antennas that formed the global VLB I array 
(VLBA-FEVN+DSN+Yg) the outer subarray of the VLA observed 
with a separate schedule, nodding between B0218-I-357 (15 mins 
on source) and 0234+285 (2 mins on source) over a period of six 
hours. The data were calibrated in AIPS using 0234-1-285 as the 
phase, amplitude and polarisation calibrator The L)-terms were 
solved for assuming a point-source for 0234+285 (a good approx- 
imation to the polarised intensity image) whilst simultaneously 
solving for the source polarisation. The position angle of polari- 
sation was again aligned using the inner VLA subarray (Yg) data. 
These data were then mapped and self-calibrated, again in AIPS. 



3 RESULTS 

3.1 Global VLBI maps (total intensity) 

The final VLBI maps are shown in Fig. ^ and are made from 
uniformly-weighted data in all cases. This results in an off-source 
rms noise of 30 /iJy beam~^ (compared to a theoretical noise of 
~10 /iJy beam~^) and a synthesised beam of 1.36 x 0.41 mas for 
the total intensity maps; these are shown in the top two panels. Here 
the images are dominated by two bright sub-components, with a 
separation of > 1 mas in both A and B. These were designated (A 
or B) 1 and 2 by Patnaik et al. (1995 ) and observations at higher fre- 
quencies (15, 22 and 43 GHz) sho w the same two sub-component 
structure (Porcas & Patnaik 1996a). The core is believed to be the 



most westerly sub-component (1) in each case, based on its com- 
pactness relative to sub-component 2. Its spectrum also turns over 
at a higher frequency, as expected for a more compact component. 
As well as the two bright sub-components though, we have also 
detected large quantities of extended emission, in both A and B. 
This, however, looks quite different in each of the images, image A 
being considerably more resolved than image B. It is exactly these 
differences that can be used as constraints on the lens model as they 
must be explicable, in the absence of effects such as scattering, by 
the action of the lensing galaxy. 

Interpreting image B is straightforward as it looks like a typi- 
cal core-jet radio source, with a knotty jet extending eastwards from 
the bright core region. This is in line with what we would expect 
given the new CJ maps of Fig. [l| However, in order to aid the eye 
in comparing images A and B, the structures that we will now de- 
scribe have been labelled on the maps in Fig. ^ The jet in B appears 
to bend, pointing in an approximate north-easterly direction imme- 
diately after emerging from the core (Jet 1), before gently curv- 
ing over the next 5 mas to point in a more south-easterly direction 
(Jet 2). Finally the jet changes direction again, returning to a north- 
easterly trajectory before fading completely (Jet 3). Returning to 
the core region, there is a small counter-component that appears to 
emerge from the south-west of B 1 . 



In image A sub-components I and 2 are, as already stated, 
clearly detected. However, they appear to sit in a large structure 
that is edge-brightened on its western side and extended to the north 
and south. This has been observed before at thi s frequency ( Kem- 



ball et al. [2001[ ) as well as at 15 GHz ( [Patnaik et al.||1995| ), but the 



greater sensitivity of our maps makes the effect much more pro 
nounced. The majority of this structure, as has been noted by pre- 
vious authors, is in fact emission from the core region that has been 
tangentially stretched along a position angle ~ —30° (perpendicu- 
lar to the A-lens galaxy separation) as part of the lensing process. 
Lens models predict just such a stretching of the brighter image in 
two-image lens systems, but it is not usually as resolved as that seen 
in B0218+357. We can quantify the magnitude of the stretching by 
noticing that the separations A1-A2 and B1-B2 are approximately 
equal (1.45 ± 0.02 mas). From observations with MERLIN and the 
VLA the flux density ratio between A and B at centimetre wave 



lengths is known t o lie between three and four (Patnaik et al 



Biggs et al 



1995; 2ohen et al 



1993 



200C). As this flux density ratio is 



due to an increase in the size of image A relative to B then image A 
must be stretched by a factor 3^ relative to B in the obse r ved d i- 
rection. A similar conclusion was reached by Patnaik et al (1995). 



The other features identified in image B are also clearly vis- 
ible in image A and with properties consistent with gravitational 
lensing. The jet emerges from the core region in A and is parity 
reversed as expected - whilst Jet 1 in B initially lies above a line 
joining Bl and B2, in A it lies below a line joining Al and A2. 
The tangential stretching displaces this jet feature away from the 
1-2 separation vector compared to B. Proceeding eastwards the jet 
emission remains tangentially stretched (the Jet 2 region is seen in 
the CJl map (Fig. |l|) as a bright extension to the east of the core) 
before fading as in B into discrete blobs. Note, however, the oppo- 
site direction of the very faint emission at the far end of the jet (Jet 
3). In A this lies along a south-east line compared to north-east in 
B. The counter-component in A is also visible, but is located to the 
north of the A1-A2 separation, as opposed to south of the B1-B2 
separation vector. 

In order to search for emission from either the core of the lens 
galaxy or a third image, we have also made maps of the region in 
between images A and B. No evidence is found for either. If we 
assume that these components would have been detected if their 
peak surface brightness was greater than the 6-g rms noise level 
of 180 /iJy beam^^, we conclude that the third image has a flux 
density less than 180 /iJy, 0.02 per cent of the component A total 
flux density. Here we have assumed that the third image will be 
unresolved as its size should be that of A (~ 10 x 10 mas) divided 
by the relative de-magnification factor of ~5000 i.e. less than the 
size of the beam. Multi-path scattering (which we show is probably 
affecting the size of image A in Section 1.2) could though broaden 



the third image and reduce its surface brightness. If this is the case 
then the upper limit given above for the flux density of the third 
image is an under-estimate. 



3.2 Global VLBI maps (polarised intensity) 



Our polarised intensity maps (P = \/Q^¥lP , bottom row of 
Fig. 0) show no polarised jet emission outside of the core region. 
Sub-components 1 and 2 though are prominent and in image A are, 
as in the total intensity maps, tangentially stretched. In image B 
polarised inter-component emission between B 1 and B2 is visible 
and appears as a small jet-like extension from sub-component Bl, 
again supporting the hypothesis that this is the core. 
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Figure 2. Global VLBI maps of JVAS B0218-H357 at 8.4 GHz. Left: B, right: A. Top: Total intensity (Stokes /). Contours are plotted at multiples (-1, 1, 2, 
4, 8, 16, etc) of 3(T where a is the off-source rms noise in the map (30 /iJy beam^^). The restoring beam is shown in the bottom-left corner and has a FWHM 
of 1.36 X 0.41 mas at a position angle of — 7?5. Bottom: Polarised intensity plotted on the same angular scale as the figures in the top row. Contours are also 
plotted at the same levels as in the top figures, but the restoring beam here has a FWHM of 1.27 X 0.38 mas at a position angle of — 8?9. 



The top row of Fig. |3| shows polarisation E-field vectors over- 
laid on total intensity contours. The effect of the differential Fara- 
day rotation between the images can be seen as a rotation of the 
polarisation vectors in image A by ~ 50° clockwise relative to im- 
age B. In image B it is possible to identify changes in the polarisa- 
tion position angle along the jet axis, including across B2, the least 
compact sub-component. Generally though the position angle tends 
to lie approximately parallel to the jet axis. This "alignment" is en- 
tirely coincidental as both the jet axis and polarisation position an- 
gles are changed from their true values, the first by the gravitational 
effect of the lens and the second by Faraday rotation. In image A 
the distribution of the position angle is more complicated, there be- 
ing an anomalous region of polarisation that stretches northwards 
from the midway point between Al and A2. Here the polarisation 
vectors are oriented at ~ 0° as opposed to elsewhere in this image 
where the majority of the polarisation is oriented between position 
an gles of ~30-40° . Thi s anomalous region is also seen in the maps 



of 



Kemball et al 



(2001) 



The bottom row of Fig. ^ shows (in greyscale) percentage po- 
larisation maps of A and B, plotted on the same angular scale as the 
polarisation position angle figures above. The inner jet of B is par- 
ticularly prominent here as it is very highly polarised, with a peak 
of 20 per cent. Interestingly, the distribution of percentage polari- 



sation looks very different in image A where the peak polarisation 
is lower at only 12 per cent. There is also no sign of a significantly 
more highly polarised sub-component between Al and A2. 



3.3 Global VLBI maps (model fitting) 

We measured the flux densities and sizes for sub-components 1 and 
2 by fitting elliptical Gaussians to the total intensity maps using the 
MPS tasks JMFIT and IMFIT. We choose to fit to the maps, rather 
than to the (u, v) data, as the biasing effect of the extended jet emis- 
sion can then be more easily reduced. This is achieved by remov- 
ing all pixels from the images whose surface brightness is less than 
2.88 mjy beam~^. This corresponds to the sixth-lowest positive 
contour in the total intensity maps of Fig. ^, above which there is 
little evidence of jet emission in either image A or B. The possi- 
bility of contamination from the irmer jet remains and the effect 
of this, although probably small due to this being relatively weak, 
will be different in each image. The positions of the Gaussians are 
constrained to lie at the maxima of sub-components 1 and 2. 

The measured flux densities and deconvolved sub-component 
sizes are shown in Table ^. Errors have been estimated by compar- 
ing the results from IMFIT and JMFIT, but are probably underesti- 
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Figure 3. Global VLBI maps of JVAS B0218-H357 at 8.4 GHz. Left: B, right: A. Top: Close up of the total intensity maps (contours) of Fig. |with polarised 
intensity E-vectors overlaid. A vector of length 1 mas equals 5 mjy beam~^ . Bottom: Percentage polarisation greyscales plotted on the same angular scale as 
the figures in the top row. For each image the darkest greyscale coiTesponds to the peak intensity in the image. 



Table 3. Results of elliptical Gaussian fitting to the VLBI maps. Tabulated 
are flux densities (Si/), deconvolved major (a) and minor (6) axes and ori- 
entation of ellipticity (</>). 



Image 


(mJy) 


a (mas) 


b (mas) 


4>n 


Al 


270 ± 5 


1.3 ±0.1 


0.5 ±0.1 


-33 ± 2 


A2 


280 ± 5 


2.3 ±0.1 


1.2±0.1 


-30 ±2 


Bl 


114 ± 3 


0.30 ± 0.02 


0.16 ±0.02 


83 ± 2 


B2 


68 ±3 


0.80 ± 0.02 


0.40 ± 0.02 


81 ± 2 



mates of the true uncertainties due to the systematic effects alluded 
to in the previous paragraph. If we consider only image B, where 



the errors are smaller, we see that sub-component 1 is the more 
compact and that its flux density i s significantly lower than that 



measured by Kemball et al. (2001). As that of B2 has remained 



unchanged over the intervening five years, these maps support the 
identification of sub-component 1 as the core and the origin of the 
variability seen in the VLA monitoring campaigns. 



3.4 VLA maps (outer subarray) 

The total intensity map of the VLA-only data is shown in Fig. ^ 
with polarisation E-field vectors overlaid. It shows all the major 
features of B0218-I-357 on arcsec scales - the radio cores A and B, 
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Figure 4. VLA (outer suban'ay) contour map of JVAS B0218-H357 at 
8.4 GHz with polarisation E-vectors overlaid. The restoring beam is shown 
in the bottom-left corner of the map and has a FWHM of 186 X 174 mas 
at a position angle of — 26?2. Contours are plotted at multiples (—1, 1, 
2, 4, 8, 16, etc) of 3(t where cr is the off-source rms noise in the map 
(50 fj,}y beam" ^ ). The inset at the top right shows the region of the Einstein 
ring without the A and B polarisation vectors. 



stein ring is much fainter. Also included in Table ^ are values for 
the total intensities, polarised intensities and polarisation position 
angles as measured from the VLBI data. These are calculated from 
the integrated I, Q and U flux densities which facilitates a com- 
parison with the values measured from the lower-resolution VLA 
(outer subarray) dataset. From the P and / measurements we cal- 
culate the corresponding percentage polarisations, m (%), which 
are also tabulated. 

The good correspondence between the polarisation position 
angles as measured from the VLA and VLBI data suggests that the 
absolute position angle calibration has been successful. Similarly, 
the excellent correspondence between the integrated polarised in- 
tensities signifies that the flux scale of the two arrays are correctly 
aligned. However, in the VLBI maps we only detect 81 per cent of 
the VLA flux density in image A and 75 per cent in image B. We 
think this is due to a combination of the poorer surface-brightness 
sensitivity of the VLBI array (most of the jet visible in Fig. |l| is 
not detected) and over-estimation of the VLA flux densities due to 
blending with the ring emission in the map. This latter hypothesis 
is supported by the fact that we detect more of the VLA flux density 
in image A (which is less blended with the ring in the VLA map) 
as well as by the very similar polarised intensities measured by the 
VLA and VLBI arrays. We have made tapered maps of the VLBI 
data in an attempt to recover more of the A and B image flux den- 
sities. More is detected (~5 per cent in each image), but the VLBI 
flux densities of A and B still fall significantly short of their VLA 
values. 



4 DISCUSSION 



Table 4. Total flux density (/), polarised intensity, (P), percentage polarisa- 
tion (m) and polarisation position angle ((/>) for images A and B as measured 
from the VLBI and VLA datasets. The P and (p VLBI values are derived 
from measurements of the integrated Q and U flux densities. The VLA val- 
ues have been calculated by model-fitting delta components to the (u, v) 
data. 



Image 


/ (mJy) 


P (mJy) 


m (%) 


</>( 


A (VLBI) 


602 


35.8 


5.9 


29 


A (VLA) 


744 


34.1 


4.6 


33 


B (VLBI) 


191 


17.1 


9.0 


83 


B (VLA) 


253 


16.8 


6.6 


82 



the Einstein ring and the radio jet. In addition, the region around the 
Einstein ring is shown with the polarisation vectors corresponding 
to images A and B omitted. This enables the polarisation vectors in 
the ring to be seen more clearly. It is interesting to see that around 
the majority of the ring the polarisation vectors are aligned radi- 
ally outwards. As gravit ational lensing preser ves the orientation of 
polarisation on the sky ( |3yer & Shavei 1992) this is somewhat un- 
expected. Although it is possible that the polarisation is Faraday 
rotated in the jet as well as in the cores, it would be surprising if 
the spatial distribution of Faraday depths was such that the polarisa- 
tion became aligned with the total intensity structure of the Einstein 
ring. However, this seems to be the case. 

We have estimated the flux densities of images A an d B by fit- 
ting t wo delta components to the {u, v) data in DIFMAP (Shepherd 



1997). The results are shown in Table ^. In fitting to the Stokes I 
data we include a third component (a ring, component type 4 of 
the DIFMAP modelfit function) to approximate the bright Einstein 
ring structure and improve the fit. No such additional component 
was added when fitting to the Stokes Q and U maps as the Ein- 



Our global VLBI 8.4-GHz maps have revealed a great deal of sub- 
structure in images A and B of JVAS B0218-I-357. Neither the lens 
galaxy core nor the third image were detected, both of which would 
be very useful constraints on the lens model, especially in this sys- 
tem where the centre of the lensing galaxy is currently uncertain.^ 
In this section we discuss three main topics: the constraints that 
the new maps can already put on lens mass models, the evidence 
that the image suffers from multi-path scattering and the different 
polarisation distributions in each image. 



4.1 Mass model constraints 

The extensive substructure that we have detected will be useful 
for lens modelling and ultimately it is our intention to use the 
LensClean algorithm on these data for this purpose, the results 
of which will be presented in a future paper. In the interim we 
have back-projected each image into the source plane in order to 
qualitatively demonstrate the accuracy of the model. This has been 
done using the CLEAN components and an isothermal model of the 
lens galaxy o ptimized using the p ositions of the 15GHz VLBI sub- 
components (Patnaik et al. 1995), a flux density ratio of 3.75 and 



the galaxy position found from our LensClean of VLA and MER- 
LIN data ( Wucknitz 2001). With x and y pointing along the major 
and minor axis of the ellipsoid, the potential is written as 



— ceo 



+ 



yZ 



(1) 



^ The lack of detection of a third image can be used though to place an 
upper limit on a finite core radi us of the lens galaxy. This has been done for 
B0218-H357 by |Norbury| ( pO(q ). 
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Figure 5. Maps of images A (right) and B (left) after being baclc-projected into the source plane. The restoring beam is shown in the bottom-left comer and 
hasaFWHMof2.61 X 0.49 mas at a position angle of-13?7. 



The critical radius is qq = 160.6 mas, the ellipticity e = 0.0583 
and the position angle of the major axis 6 — — 41.9deg, mea- 
sured from north through east. The centre of the lens is at zq — 
(260, 117.5) mas relative to A. 

In order to compensate for the very different effective resolu- 
tion of A compared to B, the back-projected CLEAN components 
of both images have been restored with the same beam. This beam 
is just large enough that both the A and B back-projected beams are 
encompassed by it, which in practice means that it is approximately 
equal to the restoring beam in the image plane back-projected from 
image B to the source plane. This makes the nominal resolution in 
each the same. The back-projected maps are shown in Fig. ^. 

Compared to the lens-plane images of Fig. ^ the source-plane 
images look quite similar. Worthy of note is the excellent corre- 
spondence between the counter-components to the west of the core 
region which are probably images of a weak counterjet. This is in- 
triguing though as the le nsed radio source in B02 18-1-357 has been 
classified as a BL Lac (O'Dea et al. 1992; ^rowne et al 



Stickel & Kuhi 



1993 



1993), a class of objects that are not expected to 



display counterjets according to the unification schemes of extra- 
galactic radio sources. This theoretical prediction is supported by 
the non-detection of counterje ts in BL Lac sources, a no table ex- 



ception being PKS 1413-1-135 ( Perlman et al. 



1994 



1996). 



A closer comparison of the A and B source-plane maps shows 
that the B jet seems to be stretched by about 10 per cent relative 
to the A jet i.e. the jet in B is longer. This stretching can only be 
explained with a different radial mass profile, compared to that used 
in the model, since the jet is directed more or less radially relative 
to the galaxy's centre. Writing the radial surface mass density as 
a function of radius as a power law, E(r) cx r"'', we find P ~ 
1.04. A very similar value is found (/9 = 1.06 =b 0.03) if we use 
the 15-GHz sub-component positions of 



constraints (Wucknitz 



20011 



Patnaik et al 



( 199^ 



4.2 Multi-path scattering 

In Fig. Is] image A looks smoother than image B. We suggest 
that this is due to scatter-broadening, which may also explain the 
frequency-dependent sizes of the images at low (^2.3 GHz) fre- 
quencies (Biggs et al., in preparation). Other examples of where 



scattering is believed to modify the surface brightness of graV' 
itationally lensed i mages are PKS 1830-211 diones et al| |l996| 



Guirado et al 



1999) and CLASS B 1933-^503 (Marlow et al. 1999) 



In B02 18-1-357, as well as in the other two systems, the scattering 
would most likely originate in the lensing galaxy as we know from 
the observed Faraday rotation and differential depolarisation that 
this is both highly ionised and non-homogeneous. Also, at the lens 
redshift of 2; = 0.6847 it is perfectly reasonable to expect that the 
ISM in front of each image, with a separation of 334 mas (2.4 kpc)^ 
between them, could be significantly different. 

In order to quantify the effect, we have fitted Gaussians to sub- 
components 1 and 2 of the back-projected images of A and B using 
the task JMFIT. The deconvolved major (a) and minor (b) axes, as 
well as the position angle of the major axis (0), are shown in Table ^ 
where it can be seen that for each sub-component the geometric 
mean of the axes (equivalent circular size) is bigger in A than in B; 
the area of the deconvolved Al is approximately eight times that 
of B 1 . Also, whilst the position angles of B are as expected (they 
point more or less down the jet as in the image plane) those of A 
remain approximately perpendicular to this. 

The advantage of model-fitting to the core sub-components in 
the source plane is that the back-projection removes the differences 
between the two images, reducing the need for removal of confus- 
ing pixels corresponding to the jet emission and rendering any bias 
that might exist identical in each. The disadvantage of fitting in the 
source plane is the potential for errors to accumulate during the 
many convolutions and deconvolutions that have been applied to 
the data at various stages of the data reduction. However, we note 
that the measured sizes in the image plane also show that image A 
is larger than image B, by a factor greater than the relative mag- 
nification. From the deconvolved sizes of Table ^ the image-plane 
area ratios Al/Bl and A2/B2 are approximately 14 and 9. For the 
remainder of this analysis we only consider the source-plane size 
measurements. 

We assume that the observed size of an A sub-component is 
the sum, in quadrature, of the intrinsic size and a scattering scale 
size {9sc)- We further assume that, for each sub-component, the true 



Throughout this paper we assume a flat universe with Hq 
70kms-i Mpc-\ Uq = 0.3 and Aq = 0.7. 



© 0000 RAS, MNRAS 000, 000-000 



Global 8.4-GHz VLBI observations ofJVAS B0218+357 9 



Table 5. Deconvolved major (a) and minor (b) axes and orientation of el- 
lipticity (<j>) of Gaussians fitted to the back-projected maps as well as their 
geometric means and the size of the derived scattering discs (9bc)- All quan- 
tities are measured in mas apart from position angles which are measured in 
degrees. We estimate eiTors on the major and minor axes of about 5 per cent. 



Image 


a 


b 


<t> 


y'a X b 


dsc 


Al 
Bl 


1.07 
0.34 


0.67 
0.28 


151 
57 


0.85 
0.31 


0.79 


A2 
B2 


1.24 
0.78 


0.96 
0.64 


149 

76 


1.09 
0.71 


0.83 



intrinsic size is actually that of its image B counterpart i.e. that this 
image is unaffected by scattering at this frequency. This assump- 
tion is supported by the similarity of the deconvolv ed size of B 1 



as mea sured from these data and the 15-GHz data of Patnaik et al , 
(1995). The derived size of ^sc is also shown in Table ^ and is the 



same for both sub-components, 0.8 mas. The true size of the scat 
tering disc (we have measured sub-component sizes in the source 
plane) is 

^sc = dsc-pr—, (2) 

where is the angular diameter distance between the Earth and 
the lensed source and Dds that between the lensing galaxy and the 
lensed source. The relation in equation ^ is analogous to that be- 
tween the tru e deflection angle (a) an d the apparent deflection an- 
gle (a) (e.g. [Refsdal & Surdej||l994| ). No correction for the lens 
magnification is required as the displacement in the image plane 
due to scattering is the scattering angle magnified by the lens, and 
this was removed when back-projecting each image. We calculate 
^sc = 3.4 mas. 

The size of the scattering disk is equal to the refractive 
l ength scale ir\ ci) subtended at the distan ce of the lensin g galaxy 
( Narayaij 1992). Using the formulation of Walkei (1998) we have 
estimated the scattering strength (SM) of the ionised medium. With 
the observed frequency of 8.4 GHz corrected to the redshift of 
the lensing galaxy we measure SM ~ 150 kpc m~^'^^^ . This is a 
very large value compared to those seen along typical lines of sight 
through the Galaxy, although comparable (and in some cases much 
higher) val ues have been measure d, particularly towards the Galac- 
tic centre ( [Cordes & Lazioj |2002| ). Therefore, although the scatter- 
ing in front of B02 18-1-357 image A is extreme, it is not ruled out 
by observations in the Galaxy. 



4.3 Depolarisation 

The polarisation properties of this system continue to intrigue. 
From observations of this system with MERLIN and the VLA, a 
general picture has emerged of image A being less polarised than 
image B and of both depolarising with increasing wavelength, al- 
beit much more steeply in A. The reason for this is probably in- 
homogeneities in the magneto-ionic medium that is responsible for 
the Faraday rotation of the polarisation position angles, inhomo- 
geneities that are greater in the region of A than B. Complica- 
tions that hinder the interpretation of the polarisation results include 
"beating" of the polarisation position angles of sub-components 1 
and 2 (for observations which do not resolve the sub-components), 
changing source structure with frequency and different magnifica- 
tion gradients across each image. 

The resolution afforded by VLBI observations allows the po- 
larisation structure of the source to be resolved and constraints to 



be put on the angular scale of the inhomogeneities in the Faraday 
screen. In Fig. ^ we see that image B, the least resolved of the im- 
ages, is seen to have a higher peak polarisation (20 per cent) than 
that of A (12 per cent). This shows that we are not seeing the effects 
of intrinsic changes in the polarisation structures across the radio 
source as the lower effective resolution of B would cause these to be 
averaged incoherently more than in A i.e. we would expect image B 
to be less polarised than image A. Therefore, our global VLBI maps 
support the theory that the depolarisation is due to non-uniformities 
in the Faraday screen which are greater in front of image A. They 
also allow us to constrain the size of these irregularities which must 
be smaller than the synthesised beam i.e. < 1 mas. This corresponds 
to a linear distance of <7 pc at the redshift of the lensing galaxy. 

We also note that the scatter-broadening identified in Sec- 
tion 12 could cause the depolarisation. This is because scatter- 
broadening effectively blurs a source, causing regions with differ- 
ent polarisations to overlap and be averaged together. As the scat- 
tering is different in front of each image the polarisations of each 
will be different and, as scattering increases at longer wavelengths, 
the observed frequency dependence of the depolarisation could also 
result. 

Another interpretation is that the differences between the im- 
ages could result from the combination of i ntrinsic polarisation 



variability and ti me delay. VLA monitoring (Corbett et al. 



1996 



Biggs et al, 1999) has shown that this source is variable in polarisa- 
tion as well as in total flux density. Two-epoch VLBI observations 
separated by the time delay, 10.5 ± 0.4 d, could help to clarify 
the polarisation structure in this source. At the same time the pic- 
ture of A being more depolarised than B is too simplistic . From 
earlier VLA monitoring observations ( [Corbett et al, 1996) it was 
found that, with the time delay removed, the polarisation of image 
B was lower than that of image A at 15 GHz, the opposite to what 
is usually observed. This could signify time variability within the 
magneto-ionic medium or movement of the source relative to this. 



5 SUMMARY AND FUTURE WORK 

We have presented new maps of the lens system JVAS B02 18-1-357 
made with a global VLBI array that combine excellent sensitiv- 
ity and resolution to produce the best maps yet of the mas-scale 
structure of this system. Unfortunately, we were unsuccessful in 
detecting the third image or the core of the lensing galaxy. With 
the new maps we have been able to explore several aspects of the 
lensing galaxy. From the different lengths of the jet in each image 
we have been able to constrain the radial mass profile and find that 
it is slightly steeper than isothermal {(3 ~ 1.04). We also find evi- 
dence for multi-path scattering and depolarisation due to the dispar- 
ity between the component sizes and polarisations in each image. 
These physical effects most likely result from propagation through 
the turbulent ionised ISM of the lensing galaxy and, in this way, 
B0218-I-357 is probably the best example of how gravitational lens- 
ing can act as a unique probe of the ISM of high-redshift galaxies. 

Th ese data will be used with LensClean ( Kochanek & 
Narayan 1l992| ) in order to exploit fully the resolved structure of the 
images and provide further constraints for the lens model. Whilst 
the effects of scatter-broadening might seem to limit the usefulness 
of these data for lens modelling, we note that the effects of the scat- 
tering and the lens can be separated. This is because scattering only 
effects the sizes of components whilst the lens changes both the 
component sizes and separations. Therefore, future modelling will 
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be able to not only constrain the lens model, but also to characterise 
the scattering more accurately than has been possible in this paper. 

Further observations with a Global VLBI array at 5 GHz may 
be worthwhile to try and detect, with higher dynamic range, the 
jet emission that has been revealed in the new CJ maps (Fig. |l|). 
The present and future VLBI maps will be used to constrain the 
lens model, in conjunction with MERLIN, VLA and HST data. 
HST observations will be carried out with the Advanced Camera 
for Surveys which our simulations predict will give us an accuracy 
of ~10 mas in the galaxy position relative to the lensed images. 
With this vital parameter in the lens model finally measured di- 
rectly (as opposed to the modelled position found from LensClean) 
the uncertainty in the i?o determination from this system will be 
reduced to a level of about 3 per cent. 
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